Independant inversion from T2 to T6ΒΆ
# -------- ResIPy log file -------
import warnings
warnings.filterwarnings('ignore')
from resipy import Project
import pyvista as pv
pv.set_jupyter_backend('panel')
API path = /home/ben/miniconda3/envs/py39/lib/python3.9/site-packages/resipy
ResIPy version = 3.3.3
cR2.exe found and up to date.
R3t.exe found and up to date.
cR3t.exe found and up to date.
pl = pv.Plotter(notebook=True,shape=(2, 3))
name = 'Salorno_2_cl5_step_rec10'
for i in range(5):
dir_simu = '../INV/regMode1step' + str(i+2) + '/'
k = Project(dirname=dir_simu, typ='R3t')
#k.loadProject(dir_simu + "Salorno_2_step" + str(i+2) + ".resipy")
k.loadProject(dir_simu + name + str(i+2) + ".resipy")
if i<3:
pl.subplot(0, i)
else:
pl.subplot(1, i-3)
pl.add_text("Time" + str(i), font_size=10)
print(i)
k.showResults(
index = 0,
attr='Resistivity(ohm.m)',
ax=pl,
vmin=0,
vmax=300,
color_map='jet',zlim=[-30,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.show()
/home/ben/miniconda3/envs/py39/lib/python3.9/site-packages/vtkmodules/util/numpy_support.py:74: DeprecationWarning: `np.bool` is a deprecated alias for the builtin `bool`. To silence this warning, use `bool` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.bool_` here.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
_vtk_np = {vtkConstants.VTK_BIT:numpy.bool,
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step2
clearing dirname
8993/9677 reciprocal measurements found.
15 measurements error > 20 %
0
/home/ben/miniconda3/envs/py39/lib/python3.9/site-packages/vtkmodules/util/numpy_support.py:74: DeprecationWarning: `np.bool` is a deprecated alias for the builtin `bool`. To silence this warning, use `bool` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.bool_` here.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
_vtk_np = {vtkConstants.VTK_BIT:numpy.bool,
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step3
clearing dirname
7998/8691 reciprocal measurements found.
17 measurements error > 20 %
1
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step4
clearing dirname
8314/9011 reciprocal measurements found.
10 measurements error > 20 %
2
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step5
clearing dirname
8329/9033 reciprocal measurements found.
16 measurements error > 20 %
3
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step6
clearing dirname
8662/9344 reciprocal measurements found.
30 measurements error > 20 %
4
for i in range(5):
dir_simu = '../INV/regMode1step' + str(i+2) + '/'
k = Project(dirname=dir_simu, typ='R3t')
#k.loadProject(dir_simu + "Salorno_2_step" + str(i+2) + ".resipy")
k.loadProject(dir_simu + name + str(i+2) + ".resipy")
pl = pv.Plotter(notebook=True)
k.showResults(
index = 0,
attr='Resistivity(ohm.m)',
ax=pl, pvslices=([5,15,25],[15],[]),
vmin=0,
vmax=300,
color_map='jet',zlim=[-30,0],
background_color='white',
pvgrid = True,
pvshow=True,
)
pl.show_bounds(all_edges=True)
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step2
clearing dirname
8993/9677 reciprocal measurements found.
15 measurements error > 20 %
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step3
clearing dirname
7998/8691 reciprocal measurements found.
17 measurements error > 20 %
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step4
clearing dirname
8314/9011 reciprocal measurements found.
10 measurements error > 20 %
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step5
clearing dirname
8329/9033 reciprocal measurements found.
16 measurements error > 20 %
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step6
clearing dirname
8662/9344 reciprocal measurements found.
30 measurements error > 20 %
pl = pv.Plotter(notebook=True,shape=(2, 3))
for i in range(5):
dir_simu = '../INV/regMode1step' + str(i+2) + '/'
k = Project(dirname=dir_simu, typ='R3t')
#k.loadProject(dir_simu + "Salorno_2_step" + str(i+2) + ".resipy")
k.loadProject(dir_simu + name + str(i+2) + ".resipy")
if i<3:
pl.subplot(0, i)
else:
pl.subplot(1, i-3)
k.showResults(
index = 0,
attr='Resistivity(ohm.m)',
ax=pl, pvslices=([5,15,25],[15],[]),
vmin=0,
vmax=300,
color_map='jet',zlim=[-30,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.show()
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step2
clearing dirname
8993/9677 reciprocal measurements found.
15 measurements error > 20 %
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step3
clearing dirname
7998/8691 reciprocal measurements found.
17 measurements error > 20 %
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step4
clearing dirname
8314/9011 reciprocal measurements found.
10 measurements error > 20 %
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step5
clearing dirname
8329/9033 reciprocal measurements found.
16 measurements error > 20 %
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/regMode1step6
clearing dirname
8662/9344 reciprocal measurements found.
30 measurements error > 20 %
Time lapse between November 2021 and April 2022ΒΆ
# -------- ResIPy log file -------
from resipy import Project
import pyvista as pv
dir_simu = '../INV/timeLapse/TL12regMode2/'
k = Project(dirname=dir_simu, typ='R3t')
#k.loadProject(dir_simu + "Salorno_TL12_reg2.resipy")
#k.loadProject(dir_simu + "Salorno_TL12_filtAppRes_estErr.resipy")
k.loadProject(dir_simu + "Salorno_TL12regMode2rec10finemeshparallelF.resipy")
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/timeLapse/TL12regMode2
clearing dirname
8993/9677 reciprocal measurements found.
15 measurements error > 20 %
7998/8691 reciprocal measurements found.
17 measurements error > 20 %
8993/9677 reciprocal measurements found.
15 measurements error > 20 %
pl = pv.Plotter(notebook=True,shape=(1, 3))
pvslices=([5,15,25],[15],[])
pl.subplot(0, 0)
k.showResults(
index = 0,
attr='Resistivity(ohm.m)',
ax=pl,
vmin=0,
vmax=300,
color_map='jet',zlim=[-30,0],
background_color='white',pvslices=pvslices,
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.subplot(0, 1)
k.showResults(
index = 1,
attr='Resistivity(ohm.m)',
ax=pl,
vmin=0,
vmax=300,
color_map='jet',zlim=[-30,0],
background_color='white',pvslices=pvslices,
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.subplot(0, 2)
k.showResults(
index = 1,
attr='difference(percent)',
ax=pl, pvslices=([5,15,25],[15],[]),
vmin=0, vmax=100,
color_map='bwr',zlim=[-30,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.show()
pl = pv.Plotter(notebook=True)
pvslices=([5,15,25],[15],[])
k.showResults(
index = 1,
attr='difference(percent)',
ax=pl, pvslices=([5,15,25],[15],[]),
vmin=0, vmax=100,
color_map='bwr',zlim=[-30,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.show()
pl = pv.Plotter(notebook=True)
k.showResults(
index = 1,
attr='difference(percent)',
ax=pl, pvcontour=[50, 100],
vmin=50, vmax=150,
color_map='jet',zlim=[-31,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.show()
pl = pv.Plotter(notebook=True)
k.showResults(
index = 1,
attr='difference(percent)',
ax=pl, pvcontour=[100, 150],
vmin=50, vmax=150,
color_map='jet',zlim=[-31,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.show()
Time lapse during tracer injection (background inv)ΒΆ
# -------- ResIPy log file -------
from resipy import Project
import pyvista as pv
dir_simu = '../INV/timeLapse/TL2b3456regMode2/'
k_back = Project(dirname=dir_simu, typ='R3t')
#k.loadProject(dir_simu + "Salorno_TL2b3456regMode2.resipy")
#k_back.loadProject(dir_simu + "Salorno_TL2b3456regMode2trimRecfinemeshparallelTlinErr.resipy")
#k_back.loadProject(dir_simu + "Salorno_TL2b3456regMode2trimRecfinemeshparallelF.resipy")
k_back.loadProject(dir_simu + "Salorno_TL2b3456regMode2trimRecfinemeshparallelF_3.resipy")
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/timeLapse/TL2b3456regMode2
clearing dirname
1072/3289 reciprocal measurements found.
2 measurements error > 20 %
1072/3289 reciprocal measurements found.
0 measurements error > 20 %
1072/3289 reciprocal measurements found.
2 measurements error > 20 %
1072/3289 reciprocal measurements found.
0 measurements error > 20 %
1072/3289 reciprocal measurements found.
2 measurements error > 20 %
1072/3289 reciprocal measurements found.
2 measurements error > 20 %
pl = pv.Plotter(notebook=True,shape=(2, 3))
for i in range(5):
if i<3:
pl.subplot(0, i)
else:
pl.subplot(1, i-3)
pl.add_text("Time" + str(i), font_size=10)
k_back.showResults(
index = 0,
attr='Resistivity(ohm.m)',
ax=pl,
vmin=0,
vmax=300,
color_map='jet',zlim=[-30,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.show()
pl = pv.Plotter(notebook=True,shape=(2, 4))
for i in range(4):
pvcontour=[0, 50]
pl.subplot(0, i)
k_back.showResults(
index = i,
attr='difference(percent)',
ax=pl, pvcontour=pvcontour,
vmin=0, vmax=150,
color_map='bwr',zlim=[-31,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.subplot(1, i)
pvcontour=[50, 150]
k_back.showResults(
index = i,
attr='difference(percent)',
ax=pl, pvcontour=pvcontour,
vmin=0, vmax=150,
color_map='bwr',zlim=[-31,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.show()
Time lapse during tracer injection (difference inv)ΒΆ
# -------- ResIPy log file -------
from resipy import Project
import pyvista as pv
dir_simu = '../INV/timeLapse/TL2b3456regMode1/'
k_diff = Project(dirname=dir_simu, typ='R3t')
#k.loadProject(dir_simu + "Salorno_TL2b3456.resipy")
k_diff.loadProject(dir_simu + "Salorno_TL2b3456regMode1trimRecfinemeshparallelF.resipy")
Working directory is: /home/ben/Documents/GitHub/BenjMy/Salorno_river/INV/timeLapse/TL2b3456regMode1
clearing dirname
1072/3289 reciprocal measurements found.
2 measurements error > 20 %
1072/3289 reciprocal measurements found.
0 measurements error > 20 %
1072/3289 reciprocal measurements found.
2 measurements error > 20 %
1072/3289 reciprocal measurements found.
0 measurements error > 20 %
1072/3289 reciprocal measurements found.
2 measurements error > 20 %
1072/3289 reciprocal measurements found.
2 measurements error > 20 %
pl = pv.Plotter(notebook=True,shape=(2, 3))
for i in range(5):
if i<3:
pl.subplot(0, i)
else:
pl.subplot(1, i-3)
k_diff.showResults(
index = 0,
attr='Resistivity(ohm.m)',
ax=pl, pvslices=([5,15,25],[15],[]),
vmin=0,
vmax=300,
color_map='jet',zlim=[-30,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.show()
pl = pv.Plotter(notebook=True,shape=(2, 2))
#pl.add_text("Diff\n", font_size=30)
for i in range(4):
if i<2:
pl.subplot(0, i)
else:
pl.subplot(1, i-2)
k_diff.showResults(
index = i,
attr='difference(percent)',
ax=pl, pvslices=([5,15,25],[15],[]),
vmin=0, vmax=10,
color_map='bwr',zlim=[-30,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.show()
pl = pv.Plotter(notebook=True,shape=(2, 4))
for i in range(4):
pvcontour=[0, 5]
pl.subplot(0, i)
k_diff.showResults(
index = i,
attr='difference(percent)',
ax=pl, pvcontour=pvcontour,
vmin=0, vmax=150,
color_map='bwr',zlim=[-31,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.subplot(1, i)
pvcontour=[5, 10]
k_diff.showResults(
index = i,
attr='difference(percent)',
ax=pl, pvcontour=pvcontour,
vmin=0, vmax=150,
color_map='bwr',zlim=[-31,0],
background_color='white',
pvgrid = True,
pvshow=False,
)
pl.show_bounds(all_edges=True)
pl.show()